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Abstract 



C^ , We present an analysis of the Feynman path centroid density that pro- 

a. 

vides new insight into the correspondence between the path integral and the 



Schrodinger formulations of statistical mechanics. The path centroid density 
is a central concept for several approximations (centroid molecular dynamics, 



^^ I quantum transition state theory, and pure quantum self-consistent harmonic 

00 

approximation) that are used in path integral studies of thermodynamic and 

cn 

o 

<^ . 

Q"^ , the quasi-static response of the equilibrium system to an external force. The 



dynamical properties of quantum particles. The centroid density is related to 



C^ I path centroid dispersion is the canonical correlation of the position operator, 

. , that measures the linear change in the mean position of a quantum particle 

Q I upon the application of a constant external force. At low temperatures, this 
O ' 

quantity provides an approximation to the excitation energy of the quantum 

•i— I , 

rS ' system. In the zero temperature limit, the particle's probability density ob- 



tained by fixed centroid path integrals corresponds to the probability density 
of minimum energy wave packets, whose average energy define the Feynman 
effective classical potential. 



I. INTRODUCTION 

The therinodynamic properties of many quantum systems have been studied within the 
path integral formulation of statistical mechanics.El An essential aspect of this formulation 
is the mapping of the quantum system onto a classical model, whose equilibrium properties 
can be derived by classical molecular dynamics or Monte Carlo simulation techniques with 
arbitrarily accuracy. However, the efficient calculation of equilibrium dynamical properties of 
quantum many-body systems by this formulation, i.e., time dependent correlation functions, 
still remains as a challenging unsolved problem.H 

Several applications of the path integral formulation have been developed around a quan- 
tity introduced by Feynman and Hibbs:El, the effective classical potential (ECP). This concept 
allows the formulation of a variational principle for quantum systems in thermodynamic equi- 
librium, whose importance can be appreciated by quoting the conclusions of the Feynman 
and Hibbs' book: " [This variational principle is] the only example of a result obtained with 
path integrals which cannot be obtained in simple manner by more conventional methods" .□ 
The original formulation of this principle used a free particle as a reference system, a fact that 
limits the range of validity of this approximation to temperatures where the system behaves 
nearly in a classical way. An essential improvement of this variational theory was formulated 
by Giachetti and TognettiB, and, independently, by Feynman and KleinertQ, using a har- 
monic oscillator as reference system. With this improvement the equilibrium properties of 
quantum anharmonic systems can be realistically approximated even in the low temperature 
limit.l3 The name "pure quantum self-consistent harmonic approximation" has been recently 
coined for the application of the variational principle at this level of approximation.^ 

Further applications of the ECP concept are related to kinetic and dynamical properties 
of quantum systems in thermodynamic equilibrium. Gillan formulated the basis of the 
quantum transition state theory (QTST), a kinetic approximation to calculate rate constants 
of thermal activated processes, that may be applied to anharmonic many-body systems at 
temperatures where quantum effects are important .&t3 This approximation has been used in 



condensed matter and chemical physics investigations, e.g. in recent studies of reorientation 
rates of hydrogen around a boron atom in doped sihconjlil, or proton transfer reactionscj. 
The most important quantity in QTST is the so called centroid density, a function that 
carries exactly the same information as the ECP: the centroid density is an exponential 
function of the ECP. An interesting dynamical approximation, called centroid molecular 
dynamics (CMD), was formulated by Cao and Vothtia to calculate real time correlation 
functions of quantum particles at finite temperatures. The ECP is an essential quantity for 
this approximation as the dynamical properties are derived from trajectories generated by 
classical equations of motion of particles moving in the ECP. A justification of the CMD 
approach has been given by showing that the centroid position correlation function agrees to 
second order, in a Taylor expansion in time, with the Kubo transformed position correlation 
function.EJ The latter correlation function is related to the dynamical response of the system 
to an external force. Unfortunately, the limits and capability of CMD for many-body systems 
at temperatures where quantum effects and anharmonicity are relevant have not been at 
present fully established. For recent applications of this dynamical approach see Refs.li^ 
andEj. 

The centroid density or, equivalently, the ECP, are then important concepts in the theory 
of path integrals. The motivation of the present work is based on a simple idea: if the centroid 
density is the basis of several interesting physical applications within the path integral 
formulation (i.e., a variational principle, and the QTST and CMD approximations), then 
the definition of this quantity within the Schrodinger formulation may lead to new physical 
insight into these approximations. Our goal is to present the correspondence between path 
integral concepts, such as the centroid density or the ECP, and the Schrodinger formulation. 
In a recent contributionli^ we have analyzed the centroid density showing that is related to 
the quasi-static response of the quantum system to a constant external force. Some physical 
implications of this analysis have been explored in the zero temperature limit, e.g., the ECP 
corresponds to the mean energy of minimum energy wave packets (MEWP, to be defined 
below) and CMD is an approximate dynamics based on these wave packets. In this work we 



present a full account of these findings, with focus on static equilibrium properties. 

This paper is organized as follows. The theoretical part is presented in Sec. |1|, which 
is divided into several subsections. Sections |11A| and [11 B| review the Schrodinger and path 



integral formulation of the equilibrium density matrix. The influence of a constant external 
force on the Euclidean action of the system is treated in Sec. [11 C| . Important path integral 
quantities, as the centroid density and the static-force response (SFR) density matrix, are 
introduced in Sec. [II D[ . In Sec. |IIE| the centroid density and the SFR density matrix are 



derived in the Schrodinger formulation by means of their moment generating functions. Sec. 



Ill presents some relevant physical consequences of the theoretical part. It includes two 



subsections. Sec. [1IIA[ deals with the static isothermal susceptibility at finite temperatures. 



and Sec [II B treats the zero temperature limit of the ECP and of the SFR density matrix. 



The concluding remarks are given in Sec. IV. 



II. THEORY 

A. Schrodinger formulation of the density matrix 

We consider the simple case of a quantum particle of mass m having bound states in a 
one-dimensional potential V{x). The extension to the many-particle case, for distinguishable 
particles, is presented in Appendix ^. For the subsequent analysis of the Feynman path 
centroid density, we study the static response of the particle to a constant external force, /. 
Therefore, we consider the Hamiltonian of the particle, H{f), as a function of the external 
force: 

H{f) = H-fx, (1) 

where x is the position operator and H = H{0) is the Hamiltonian operator in the 
absence of the force: 

6 = ^^V(^), (2) 



with p being the inomentuin operator. The following convention is used throughout the 
text: operators (or functions) that share the same name, as H{f) and H , but differ either 
by the presence or number of arguments, represent different operators (or functions). The 
eigenfunctions, |^„(/)), of the Hamiltonian, H{f), satisfy the Schrodinger equation: 

H{f) \Mf)) = EnU) \Mf)) , (3) 

where En{f) are the corresponding eigenvalues. In the position representation, the eigen- 
functions are expressed as: 

Mx; f) = {x\Mf)) , with n = 0, 1, 2, ... . (4) 

If we consider a canonical ensemble of independent particles in thermal equilibrium at 
temperature T, the statistical state of the system in the presence of the external force, /, is 
described in terms of the unnormalised density operator: 

Pif) = e-'^^^^ , (5) 

where /3 is the inverse temperature {kBT)~^, with ks being the Boltzmann constant. 
The position representation of this operator is: 

oo 

p{x,x'; f) = {x\p{f)\x') = E e-'^^"(^Vn(x; /)C(x'; /) . (6) 

n=0 

The unnormalised particle's probability density is given by the diagonal elements of the 
density matrix, p{x,x; /). The canonical partition function in the presence of the external 
force is defined by the trace of the density matrix: 

/oo °° 

rfxp(x,x;/) = Ee"^''"^^^ (7) 

-~ n=0 

where the second equality is obtained by substitution of p{x, x; f) by the expression given 
in Eq. m. 



B. Path integral formulation of the density matrix 

In the path integral formulation the elements of the unnormalised density matrix in the 
presence of the external force are given byH 

px'=x(f3h) 

p^x,x';f)= D[x(u)]e-^^["(")'^1 . (8) 

Jx=x{0) 

u is an imaginary time that varies between and j3?i, and S[x{u)] /] is the functional of 
the Euclidean action of the path x{u): 

S[x{uy, f] = j^ du (^-x{uf + V[x{u)] - fxiu)j . (9) 

The function x{u) represents the derivative of x{u) with respect to u. The integral 
measure of the path integral is given by: 

zjixMi = jto n d.. (^5^] , (10) 

where the path x{u) has been discretized as (x,a;i,X2, ...,xn-i,x')B 

C. Euclidean action under an external force 

It is important to note the effect of a constant external force on the Euclidean action of 
a given path x{u). As the external force is independent of the imaginary time u, it can be 
taken out of the time integral in Eq. (P) with the result: 

S[x{uyj] = S[xiu)]-f/3hX, (11) 

where ^[^(m)] is the Euclidean action in the absence of the force, and X is the average 
point of the path, or path centroid: 

1 /•/3'i 

X = — dux{u) . (12) 

ph Jo 

For the set of paths contributing to the path integral in Eq. (^, the relation of having the 
same path centroid is an equivalence relation. Then, the sum over paths may be decomposed 
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into a sum over equivalence classes.El A class of paths is the subset formed by all paths that 
have the same centroid coordinate X. We note that if two paths, say x{u) and x'(u), belong 
to the same class, then the contribution of the external force to the Euclidean action, [i.e., 
the term —fjShX in Eq. (|11])], is identical for both paths. In other words, the contribution of 
a constant external force to the Euclidean action is identical for all paths that have the same 
centroid. This simple property is the origin of the interest of fixed centroid path integrals 
in statistical physics. 

D. Static-force response density matrix and centroid density 

The most important quantities to be defined for the class of paths with centroid at X 
are the SFR density matrixa and the centroid densityll^. We have introduced a new name, 
static-force response density matrix, because we will show that this matrix depends on the 
response of the system to constant external forces of arbitrary magnitude. The name reduced 
density matrix has been used before for this quantity.^ However, this name is less convenient 
as it is often encountered in another context. The SFR density matrix is a constrained path 
integral over a given class of paths, defined by introducing a delta function in the integrand 
ofEq. i): 



a(x,x 



rx' ( 1 /"/ST!. \ 

';X;/) = y^ B\x{u)\b\X-—\^ rf«x(«)j e-^^[-(")^^] . (13) 



The centroid density, C{X\ /), in the presence of the external force /, is defined as the 
trace of the SFR density matrix: 

/oo 
dxa{x,x-X-f). (14) 

-oo 

By introducing the expression of the Euclidean action, [Eq. ([TT|)], into Eq. (p!3D, we 
derive an equation that relates the SFR density matrix in the presence of the external force 
/, with the same quantity in absence of the force: 

a{x,x'\X-J)=o{x,x';X)e<^^^ . (15) 



By setting x = x' in the last equation and integrating over x, we get the following 
analogous relation for the centroid density: 

C{X; f) = C{X) e^f"" . (16) 

The last two equations show that if the centroid density, C{X), and the SFR density 
matrix, a{x,x';X), are known for a given class of paths in absence of an external force, 
then the corresponding quantities in the presence of the force / are easily obtained by 
multiplication with the constant factor e^^^ . Dividing the left and right hand sides of the 
last two equations, we get: 

[C{X; /)]- V(x, x'; X; /) = [C(X)]- V(x, x'; X) . (17) 

This equation implies that the normalized SFR density matrix defined for a given centroid 
position X is an invariant quantity with respect to the value of the force /. The normalization 
constant is the trace of the matrix, i.e., the centroid density. Setting x = x', we get: 

<j>'{x; X) ^ [C{X- /)]- V(x, X- X- f) . (18) 

0^(a;;X) represents, as a function of x, a normalized particle's probability density. The 
normalization of this probability density is easily checked by integrating Eq. (|l^) with 



respect to the variable x, and with the help of the definition given in Eq. (|T^ 

(ix02(x;X) = 1 . (19) 



This relation is valid irrespective of the value of the centroid coordinate, X. 

The results of this Section may help to understand path integral results that are usually 
presented from a different point of view. For example, the CMD approximationS is an ap- 
proximate dynamics for a normalized density matrix whose diagonal elements in the position 
representation are given by (j)'^{x\X). An important property of this density matrix is its 
invariance with respect to the application of an external force to the particle. Associated to 
each centroid position, X, there is a different density matrix, and the approximate dynamics 
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of CMD is formulated through a dynamical equation for the time evolution of the centroid 
coordinate. 

We are now in position to relate the SFR density matrix and the centroid density to 
physical observables, an essential step for the understanding of important path integral 
concepts from the point of view of the Schrodinger formulation. 

E. Moment generating function of the centroid density 

After the decomposition of the path integral given in Eq. (||) into disjoint classes, we can 
recover the whole path integral, which is the quantity related to physical observables, by an 
integral over all classes. The elements of the unnormalised density matrix, in the presence 
of an external force /, are then obtained as: 

/oo 
dXa{x,x';X;f). (20) 

-oo 



Introducing in this equation the expression given in Eq. ([15D for the SFR density matrix, 
one gets: 

/CO 
dX(T{x,x'-X)e^f^ . (21) 

-oo 

Setting X = x' and integrating over the variable x [with the help of Eqs. (|^) and (0)], 
one arrives at: 

/oo 
dX C{X) e^^^ . (22) 

-oo 

The last two equations provide the essential link for the correspondence between fixed 
centroid path integrals and the Schrodinger formulation. Although Eq. ( [2^ ) can be found 
in Kleinert's book,t3 the physical implications of this relation have not been explored until 
recently.ll3 The physical content of Eqs. (21) and (^) can be derived by two complementary 
points of view, as suggested by the mathematical structure of these equations: 

• The centroid density, C(X), in absence of the external force, is related to the parti- 
tion function, Z{f), by an integral transformation defined by the kernel e^-^"^. The 
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same relation holds between the SFR density matrix, a{x, x'; X), and the density ma- 
trix p{x,x'; /). This integral transformation is performed with respect to the centroid 
variable X, and the product jSf is the corresponding conjugate variable. Note that 
a general property of an integral transformation is that the physical information con- 
tained in the direct and transformed functions is identical. Therefore the centroid 
density C{X), as a function of X, carries out the same physical information as the 
partition function Z{f), as a function of /. The same relation holds for (t{x,x';X) 
and p{x,x'; /). This integral transformation is called a two-sided Laplace transform, 
with properties similar to those of a Fourier transform.EJ 

• The relation between the centroid density, C{X), and the partition function in the 
presence of a force, Z{f), can be interpreted as the relation between a probability 
density for the variable X ant its moment generating function. The same relation 
holds for a{x, x'; X) and p(x, x'; /). In the following, we develop this point of view. 

The moments of the centroid density in the presence of an external force are defined as: 

The brackets {•••}/■ indicate an average over the centroid density, C{X; /), and the second 
equality is obtained from Eq. (^). The moments of X in absence of the external force are: 

/oo 
dX C{X) X" , (24) 

-oo 

where the partition function is Z = Z{0). Note that the subindex "/" is omitted from 
the brackets when / = 0. The definition of the moment generating function|l3 M(/5/), of 
the normalized centroid density, Z~^C{X), is: 

e/3/^l = Z-^ / dXC{X) e^^^ . (25) 



From Eqs. (|22D and (^) we get: 



M{Pf) = ^ . (26) 
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This result has a clear physical meaning: the ratio of partition functions Z{f)/Z is the 
function generating the moments of the centroid density. Therefore, these moments may be 
defined as:" 



and with the help of Eq. (p6|) , we get: 



d{Pf) 



(27) 



/=o 



Z/3" 



g"^(/) 

a/- 



(28) 



/=o 



The l.h.s. of the last equation, i.e., the moments of the centroid density, is a quantity 
typically defined by fixed centroid path integrals, while the r.h.s. contains physical quan- 
tities, that are defined within the Schrodinger formulation. It is more interesting to study 
the cumulant generating functionJlaH^I K{Pf), of the centroid density, which is defined as 
the logarithm of the moment generating function: 



ir(/3/) = ln^ = -/3[F(/)-F] 



(29) 



where F{f) is the free energy derived from the partition function: Z{f) = exp[—(3F{f)] 
and F = -F(O). The cumulants, «:„, of the centroid density are obtained as: 

'dK{(3fy 



diPf) 



(30) 



J/=o 



With the help of Eq. (H) we get: 



P' 



n—l 



df" 



(31) 



J/=o 



The last equation shows that the cumulants of the centroid density are related to the 
change in the free energy as a constant external force is acting on the particle. The first 
cumulants of the centroid density correspond to the mean value and the dispersion of the 
centroid coordinate: 



«:i = {X} 



df 



(32) 



J/=o 



11 



n, = 6X^ = {X^} - {Xr = -\ 



d'F{f) 



op 



(33) 



/3 

The centroid dispersion, (5X^, has been related with a "classical delocalization" of the 
particle." The precise physical meaning of this quantity will be presented in Sec. PTT[ 

As the mathematical structure of Eqs. (2T) and ( P^j) is identical, the derivation done 
for the moment and cumulant generating functions of the centroid density can be repeated 
step by step for the SFR density matrix. The most important result of this analysis may 
be enunciated as follows: the cumulant generating function, Kcr{l3f), of the SFR density 
matrix, a{x,x';X), is the logarithm of the ratio of the elements of the density matrix in the 
presence and in absence of the external force f: 

^"W)-'»^^- (34) 

As an application of the analysis presented so far we derive in Appendix ^ the centroid 
density and the SFR density matrix for a harmonic oscillator by means of their cumulant 
generating functions. The information needed for this task is the value of the partition 
function, Z{f), and the density matrix, p{x, x'; /), as a function of the external force /. This 
example shows how results previously obtained by solving fixed centroid path integralg3 can 
be easily derived in the Schrodinger formulation. 

To close this theoretical Section, we explain the precise meaning of considering the cen- 
troid density as a classical-like density for the quantum particle, which is one of the suggestive 
pictures derived from the QTST and CMD approximations. Let p'^'"(x) be the classical limit 
of the unnormalised particle's probability density in absence of an external force. Then the 
classical partition function Z^'"(/), in the presence of a constant external force, is given by: 

/oo 
dxp''^\x)e'^f'' . (35) 

-oo 

This equation has formally the same structure as Eq. ( p2D for the centroid density, and 
illustrates in which sense the centroid density behaves as a classical density for the quantum 
particle. In the classical limit the ratio Z^'"(/)/Z'^'" is the moment generating function of 
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the classical particle's probability density, while in the quantum case, the ratio Z{f)/Z is 
the moment generating function of the centroid density. 

The unnormalised particle's probability density is given in the classical case by a function 
of the potential energy: 



■,'^l^(^\ - I ^ \ ' ^^t3V{x) 



In analogy to the classical result, the ECP, Fecp{X), for the quantum particle is defined 
from the centroid density as: 

We remark that the ECP is sometimes called the potential of mean force,!!^, the centroid 
potential,t3 the effective centroid potential,t3 or the quantum effective potentiahEll However, 
these names refer all to the same physical quantity. 

III. PHYSICAL IMPLICATIONS 

In this section we focus on some physical implications readily derived from the corre- 
spondence between fixed centroid path integrals and the Schrodinger formulation. Firstly 
we present some results valid at arbitrary temperatures. 

A. Static susceptibility 

The first moments of the centroid density have a clear physical meaning. We recall that 
the average position of the quantum particle in thermal equilibrium in the presence of an 
external force, /, is given by: 

/oo 
dxxp{x,x-J). (38) 

-00 



From the definition of the SFR density matrix and of the centroid density [Eqs. (|T3[) and 
(|I^)], it is easy to show that the average centroid position, defined from Eq. ( P^D by setting 
n = 1, and the average particle position are identical quantities: 
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s(/) = {X]f 



dF{f) 
df 



(39) 



where the last equahty is consequence of an exact generahzation of Eq. (|3^) , derived for 
/ = 0, to an arbitrary value of the external force. From Eqs. (|33|) and (|39|), we obtain the 
following expression for the dispersion of the centroid coordinate: 



sx^ = {x^} - {xy 



/3 



dxjf) 
df 



(40) 



J/=o 



The derivative that appears in the last equation is the static isothermal susceptibility, 
Xxxy ^ quantity that describes the static response of the average position of the quantum 
particle, with respect to the application of an external force. BJ'CJ We have then: 



6X' 



(3 



T 
/(-xx 



(41) 



This isothermal susceptibility, xix-! i^ related to the canonical correlation oi the position 
operator, (x; x) by the following relation:E3 



xlx = P {{£; £) - x'^) 



(42) 



where x = x{0) and the canonical correlation is defined as:c3 

f/3 



{x;x) = Z-'/3 



^'^-^ ' rfATr 





exp{—(3H) exp{XH)x exp{—XH)x 



(43) 



Comparing the Eqs. ( p|) and (^21) , and noting that x = {X}, we deduce that the second 
moment of the centroid density is identical to the canonical correlation of the position 
operator: 

{X'} = {x-x). (44) 

This result is consistent with the analysis presented in Ref.llj concerning the relation be- 
tween the centroid time correlation function and the Kubo transformed position correlation 
function. Eq. ( Pj) is interesting not only to understand the physical meaning of the centroid 
density, but also for practical applications. In the case of a harmonic oscillator of angular 
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frequency u the dispersion of the centroid coordinate, (5X^, is given in Eq. ( |B5|) . The static 
isothermal susceptibihty for the harmonic case is then: 

For an arbitrary anharmonic potential, we expect that at low temperatures the static 
isothermal susceptibility will be determined by a frequency close to the first excitation 
energy, A£^, of the system. Then, by substitution in the last equation of u by AEapp/h, we 
get the following approximation to AE: 

^^-^^- = K^)' ''" 

The capability of this approximation has been checked as a function of temperature for 
several one-dimensional model potentials, that are listed in Table |. The potentials V2, V4, 
and Vio are power functions of x, whose coefficients were chosen so that a particle of mass 
16 au displays the same value of the excitation energy in each of these potentials. Vdw is a 
double-well potential where the first excitation energy corresponds to the tunnel splitting. 
By Monte Carlo path integral simulationsES we have obtained the value of the centroid 
dispersion, SX"^, as a function of temperature, and the excitation energy has been estimated 
by Eq. (^6]). The results are shown in Fig. |1|, where the value of AEapp and the thermal 
energy, ksT, are displayed in units of the AE associated to each potential. For the harmonic 
potential, V2, the approximation is exact at all temperatures. For the potentials V4 and \4o 
the approximation is remarkably good at temperatures where the thermal energy is lower 
than about 1/4 of the first excitation energy. In this temperature range the tunnel splitting 
energy is approximated with an error of about 25 %, a value that at least is of the correct 
order of magnitude. This approximation can be applied to many-body problems, using the 
centroid dispersions resulting from the diagonalization of the tensor given in Eq. (|A3|) . 
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B. r ^ limit 

In the zero temperature limit, the ECP and the normahzed particle's probability densi- 
ties, 0^(x; X), [defined in Eq. ([18D ] have a simple physical meaning: they are related to the 
eigenvalues and eigenfunctions of the Hamiltonian H{f). 

We are going to combine two pieces of information to obtain the low temperature limit 
of 0^(x;X). Firstly, we know that, as T — > 0, the normalized particle's probability density 
derived from the path integral in Eq. (^) converges towards the probability density of the 
ground state of the Hamiltonian H{f). The average value of the position operator x for 
this ground state is xo(/). Secondly, we have derived that the dispersion of the centroid 
coordinate goes to zero in this limit [see Eq. (^O])]. This implies that the centroid density, 
C{X] /), must be a delta function of X centered at its average value {X} = Xo^f) [see Eq. 
(^)]. Therefore only the class of paths with centroid coordinate at X = Xo(/) contributes to 
the path integral in Eq. (H). As a consequence, the normalized probability density, 0^(x; X), 
derived for the class of paths with centroid at X = xo{f) should be identical to the ground 
state probability density of the particle in the presence of the external force /: 

hm02(a;;X) = \Mx;f)\' ,for X ^ xo{f) . (47) 

As illustration of this result we display in Fig. g the probability densities, 0^(x;X), 
obtained by fixed centroid Monte Carlo simulations of a particle in the model potentials V2, 
V4, and Vio at a temperature of ksT = 0.001 au, which is about 1/300 of the first excitation 
energy, and therefore a good approximation to the low temperature limit. For each model 
potential, the simulations were performed at three different centroid positions, X. The 
ground state densities |^o(a;;/)P obtained by numerical solution of the time independent 
Schrodinger equation, for different values of the external force, /, are also shown in the 
figure. Both probability densities are identical, apart from tiny deviations (not visible in 
the scale of the figure) due to the finite temperature and the statistical uncertainty of the 
simulation. In the harmonic case the displayed probability densities, 0^(a;;X), differ only 
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by a rigid displacement as a function of the centroid position X. We note that these curves 
are identical to the probability densities of the coherent states of the harmonic oscillator .Ej 
In the following we show that the functions (f?{x\ X) are related to a variational principle. 
Suppose that we look for a quantum state, |\1/), of the particle with Hamiltonian H whose 
mean energy, 

E,mn{X) = mH\^) , (48) 

is minimum against small variations of |\1/). Moreover, this state must satisfy two con- 
straints: its mean position is fixed at an arbitrary value X, and it is normalized: 

X = {^\x\^) , (49) 

1 = (^1^) . (50) 

By straightforward application of calculus of variations (see Appendix |C|) , one finds that 
|\1/) must be the ground state of a Hamiltonian H{f), i.e., |\1') = |V^o(/)): 

{H - fx)\Mf)) = Eo{f)\Mf)) , (51) 

where the force / and the corresponding ground state energy, -Eo(/); of the Hamiltonian 
H{f), appear as Lagrange multipliers. The value of / must be chosen so that the constraint 
in Eq. (^Q]) is satisfied. The fixed mean position X corresponds then to the average position 
of the ground state lipoif)), i.e. X = Xo{f). The average energy, Emin{X), of these minimum 
energy states is derived as a function of X with the help of Eqs. (0) and (^91) as: 



E™„(X) = Eoif) + fX Jot X = xo{f) . (52) 

We call the states |^o(/)) the MEWP's of the unperturbed Hamiltonian Hi^ We can now 
reinterpretate Eq. ( ^Tf ) by saying that, in the zero temperature limit, the probability density 
(f)'^{x;X) corresponds to the MEWP's of the Hamiltonian H. A characteristic property of 
these states is that their average energy is stationary (i.e., minimum) with respect to any 
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arbitrary change in their probabihty density that leaves constant their average position X. 
In the following we show that this average energy is identical to the ECP. 

As T — > we have derived that the centroid density, C{X; /), is a delta function centered 
at the average position of the ground state of H{f), i.e. {X}f = xo{f). We also know that 
the integral of C{X; f) with respect to X is identical to the partition function Z[f) [see 
Eqs. ( p2|) and ([T6|) ]. Then, in the zero temperature limit we can write: 



^mCiX;f) = Zif)5[X-Xoif)]. (53) 

The asymptotic behavior of C{X] /) as T — > is given by an exponential of the ECP 
exp[—PFecp{X; /)], while the asymptotic behavior of Z{f) is given by an exponential of the 
ground state energy exp[— /JE'ol/)]- From Eq. (0) the asymptotic behavior of C{X; f) and 
Z{f) should be the same at X = Xo{f), that is: 

hm F,,,(X; /) = Eoif) M X = Xo{f) . (54) 

From the property given in Eq. (|TB|) for the centroid density, and from the definition of 
the ECP in Eq. ( |57| ) it is easy to derive a relation valid at arbitrary temperature: 

Fe,p(X;/) = Fe,p(X)-/X. (55) 

By substitution of the last expression into Eq. (|54D we obtain the desired result: 



lim F,,p(X) = Eoif) + fX = ^™„(X) , for X = a;(/) , (56) 



where the last equality corresponds to Eq. ([5^). Our final conclusion is that, in the limit 
T — > 0, the value of the ECP, Fecp(X), is equal to the average energy of that MEWP whose 
average position coincides with the centroid coordinate X. The MEWP's turn out to be 
the ground states of the Hamiltonian H{f). The results obtained in this Subsection may be 
seen as a consequence of the original variational principle of Feynman and HibbsS, when it 
is applied to the zero temperature limit. However, the important role of the external force 
does not appear in the original path integral formulation of this variational principle. 



The CMD approximation is, in the T ^ hniit, an approximate dynamics based on 
MEWP's, which is accurate even for highly anharmonic potentials.Ej Cao and Voth have 
shown that CMD reproduces correctly the classical limit at high temperatures.li^ The classi- 
cal limit of our formulation of CMD in terms of MEWP's can be derived by the generalization 
of the results at T ^ to arbitrary temperatures using the MEWP's given by the excited 
states |\1'„(/)) (see Appendix P). From this generalization we find that CMD can be for- 
mulated for harmonic systems at arbitrary temperatures as an exact MEWP's dynamics, 
even in the classical limit .E3 The most important applications of CMD are related to con- 
densed phase quantum dynamics of anharmonic systems-EJO'Ea From our analysis in terms 
of MEWP's we find that for anharmonic potentials CMD provides accurate results in both 
the T ^ and the classical limits, but further work is needed to clarify the capability of 
CMD at intermediate temperatures. 

IV. SUMMARY 

We have presented the correspondence between fixed centroid path integrals and the 
Schrodinger formulation. This analysis shows that the Feynman path centroid density and 
the SFR density matrix depend on the static response of the quantum system to a constant 
external force. The path centroid density is related by a simple integral transformation to 
the partition function of the quantum system under the action of an external force. The same 
integral transformation is found between a classical phase space density and the classical 
partition function. Therefore, the interpretation of the centroid density as a classical-like 
density for the quantum system, which is one of the ideas suggested by the QTST and 
CMD approximations, has a precise physical meaning: the path centroid density behaves 
"classically" in the sense that under the action of an external force it has the same static 
response properties as a classical phase space density. The same integral transformation is 
found between the normalized canonical density matrix in the presence of an external force 
and the SFR density matrix. Within the Schrodinger formulation it is essential to introduce 
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a constant external force to define the centroid density and the elements of the SFR density 
matrix as transformed functions. However, in the path integral formulation the introduction 
of external forces is not needed, because one works directly with the transformed functions. 
This facility of the path integral formulation is the origin a large number of path integral 
applications, from a variational approximation of the thermodynamic properties of quantum 
system,0 to the CMD approximation of time correlation functions of quantum particles in 
thermal equilibrium.Ej 

In the present work, we have not tried to present the implications of our formulation in 
current applications based on fixed centroid path integrals. Nevertheless, our analysis has 
led to results that clarify the physical meaning of fixed centroid quantities. In particular, we 
have shown that the dispersion of the centroid coordinate is related to the static isothermal 
susceptibility, and that, in the zero temperature limit, the fixed centroid path integrals are 
related to the MEWP's of the unperturbed system. At T = 0, the normalized SFR density 
matrix is a pure state density matrix corresponding to a MEWP of the Hamiltonian, and 
the Feynman effective classical potential is the average energy of the MEWP's. 
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APPENDIX A: EXTENSION TO MANY-PARTICLES 

The results presented for a particle moving in one dimension can be easily generalized to 
multidimensional n-body systems of distinguishable particles. We denote the set of n particle 
coordinates as a 3ra-dimensional vector r = (ri, r2, ..., rn) where the i-particle position vector 
is Ti = {vix, Viy, Tiz). The set of external forces is f = (fi, f2, ..., fn)- The external force fi acts 
on the particle ri through the linear term, — fifi, appearing in the potential energy of the 
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Hainiltonian H{i). The centroid density associated to the n-body Hamiltonian in absence 
of external forces, H, can be defined by a generahzation of Eq. (l2^ ) as: 



... / dr C(X) e^*^ 



(Al) 



where X = (Xi,X2, ...,Xn) is a 3n-diniensional vector formed by the centroid positions of 
each particle. The last equation represents a 3n-diniensional two-sided Laplace transform. 
The cumulants of the centroid density are obtained as the derivatives of the free energy -F(f ) 
associated to the partition function Z{f): 

'dF{fy 



{^ix} 



dk 



Jf=0 



{XixXjy} - {Xix} {Xjy} 



/3 



dfixdfj 



jy 



(A2) 



(A3) 



f=0 



The last expression represents the components of a tensor that is related to the static 
isothermal susceptibity tensor by multiplication by the constant f3. 

APPENDIX B: CENTROID DENSITY AND STATIC-FORCE RESPONSE 
DENSITY MATRIX FOR A LINEAR HARMONIC OSCILLATOR 

We consider a canonical ensemble of particles of mass m moving in a one-dimensional 
harmonic potential Vho{x) = {l/2)muj'^x^. In absence of an external force, the partition 
function and the free energy area 



2sinh 



m 



F = kBTln 



2sinh 



'phcu' 



(Bl) 



(B2) 



By application of a constant external force, /, the potential changes to Vho{x) — fx. The 
new potential energy is also a quadratic function of x, and the corresponding free energy is 
easily derived as: 
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F{f) = F- 



P 



(B3) 



From Eqs. (^) and (|33|) we obtain the first two cumulants of tlie centroid density C{X) 



as: 



{X} 



f 



muo^ 



(B4) 



J/=o 



5X^ 



(B5) 



Tlie liiglier order cumulants, wliicli are proportional to successive derivatives of the last 
equation with respect to /, are zero. This result implies that the centroid density must be 
a Gaussian function of X: 



C{X) = NiGx{X-5X^ 



(B6) 



where Ni is a normalization constant and Gx{X\ 5X'^) is a normalized Gaussian function 
of X: 



Gx{X; 5X 



1 



2'k5X^ 



exp 



{X-Xf 
25X^ 



(B7) 



The constant A^i can be obtained from Eq. ( ]22|) by setting / = 0: 



/oo 
dXC{X) = Ni 



(B8) 



The centroid density of the harmonic oscillator is then: 

1 \ 



C{X) = ZGx{0; 



Pmuj'^ 



(B9) 



The last result has been derived several times using fixed centroid path integrals-B'ty'tZltHl 
The ECP for the harmonic oscillator is obtained from its definition in Eq. (^) and the last 
equation as: 

1 



Fecp{X) = F-kBT In {(3nuj) + ^muj^X^ 



(BIO) 
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The SFR density matrix a{x,x';X) may be derived following the same scheme. The 
cumulant generating function in this case was given in Eq. (p^. The elements of the 
density matrix for the harmonic oscillator in absence of an external force are:c3 



p{x,x') 



niuo 



2Tch smh.{f3huj) 



exp 



Tnuj , r, ,2 



2h smh{f3huj) 



{x^ + x'^) cosh{phuj) - 2xx' \ . (Bll) 



While in the presence of an external force we obtain from the definition in Eq. 



2muj'^ 



p(x, x'] f) = exp ( ^_ 2 P(^ " ^™"' ^' ~ ^™«) ' (B12) 



where we have used the following relations: 



2muj'^ 



Enif) =En-^,, (B13) 



i^n{x] f) = ^Pn{x - X^in) ■ (B14) 

Xmin = f /{'m^'^) is the positiou of minimum energy for the potential Vho{x) — fx. The 
cumulants of the centroid variable X, with a{x, x'; X) as its probability density, are evaluated 
as the derivatives of the cumulant generating function [Eq. p^)] with respect to the force 
at / = 0. We find that only the first and second cumulants (ki, k,2) are different from zero. 
Therefore a{x, x'; X) must be a Gaussian function of the variable X, with ki being the mean 
value of X, and H2 its dispersion: 

a(x, x'; X) = N^Gxi'^i; ^^2) ■ (B15) 

The result for the first cumulant is obtained from the cumulant generating function after 
straightforward algebra as: 

and the result for the second cumulant is: 

K, = ^{6x^-5x1^). (B17) 
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The constants 5x^ and fe^^^ correspond to the dispersion of the position coordinate for a 
canonical ensemble of harmonic oscillators in the quantum and classical case, respectively: 

5x^ = ^coi\.(^-^\ , (B18) 



The normalization constant, N2, is determined from Eq. ( pTD by setting / = 0. We get: 

N2 = p{x,x'). (B20) 

The final result is: 



a{x, x'; X) = p{x, x')Gx 



^^rln I X -\- X \ OX^j„ / (- 2 r 2 \ 

[OX — Qx^iJ 



(B21) 



5x2 V 2 / ' 5x2 

The elements of the SFR density matrix are the product of the elements of the canon- 
ical density matrix by a normalized Gaussian function of the centroid coordinate X. This 
equation is identical to Eq. (B.4) of Ref.a (apart from a different grouping of factors). 

APPENDIX C: MINIMUM ENERGY WAVE PACKETS 

We want to find the MEWP's associated to the Hamiltonian H. To simplify the notation 
we use the following abbreviations: ^ = (a;|\E'), \1'' = d^/dx, ^" = d'^'^/dx'^, V = V{x), 
and D = -h'^/{2m). 

The MEWP's minimize the functional: 



E™„ = / dx (D^^" + V^2) (CI) 

J —00 

with the constraints given in Eqs. ( ^ , |50D for the wave function \E'. We call: 



where / and E are two Lagrange multipliers. From the Euler equation: 
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one derives the differential equation tliat must be satisfied by tlie MEWP's: 

D^" + {V - /x)* = E^ , (C4) 

wliicli is tlie time independent Hamilton equation corresponding to the Hamiltonian 
H{f)- AH functions satisfying this equation, i.e., the eigenfunction |\1'„(/)), are MEWP's 
in the sense that they minimize the functional in Eq. (|Cl|) subject to the constrains given 
in Eqs. (|^,0). The excited states {n > 1) correspond to local minima of the functional. 
In the study of the T —>■ limit of the SFR density matrix we use only the states |\t'o(/)) 
corresponding to the global minimum of the functional, but the excited states are necessary 
to generalize this study to arbitrary temperatures. 
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FIGURES 
FIG. 1. First excitation energy estimated from the dispersion of the centroid coordinate for 

four model potentials at different temperatures. For anharmonic potentials the approximation 

works only at low temperature. For each model potential the true excitation energy, AE, was used 

as energy unit. 

FIG. 2. Normalized probability densities, (p'^{x;X), obtained by fixed centroid path integral 
simulations at fe^T = 0.001 au, using the model potentials V2, V4, and Flo- For each potential 
model, the centroid coordinate was fixed at X = —0.5, 0, and 0.5 au, respectively. Superposed to 
the simulation results there are the probability densities of the ground states of these potentials in 
the presence of an external force /. The value of / associated to each ground state is given in au. 
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TABLES 
TABLE L One-dimensional model potentials. The lowest energy eigenvalues (£"0, Ei, and E2) 

are given for a particle of mass 16 au. The first excitation energy is AE = Ei — Eq. The values 

were obtained by numerical solution of the time independent Schrodinger equation. All units in 

au. 
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AE 
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